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ABSTRACT 

We consider the astrophysical implications of the diffusion equation solution 
for Ultra-High En ergy Cosmic Rays (UHECR ) in the expanding universe, ob- 



tained in paper I (IBerezinsky fc Gazizovll2006l ). The UHECR spectra are calcu 



lated in a model with sources located in vertices of the cubic grid with a linear 
constant (source separation) d. The calculations are performed for various mag- 
netic field configurations {Be, Ic), where Ic is the basic scale of the turbulence and 
Be is the coherent magnetic field on this scale. The main purpose of these cal- 
culations is to demonstrate the validity of the solution obtained in paper I and 
to compare this solution with the Syrovatsky solution used in previous works. 
The Syrovatsky solution must be necessarily embedded in the static cosmolog- 
ical model. The formal comparison of the two solutions with all parameters 
being fixed identically reveals the appreciable discrepancies between two spectra. 
These discrepancies are less if in both models the different sets of the best-fit 
parameters are used. 

Subject headings: extragalactic cosmic rays, diffusive propagation of cosmic rays. 



1. Introduction 



Diffusive propagation o f ultra-high energy co s mic rays (UHECR) in extragalactic space 
been recently studied by lAloisio fc BerezinskyI (120041 . l2005l ): lLemoind (120051 ) : lAloisio et al. 



- 2 - 



( 2007a ) using the Syrovatsky solution (see Syrovatskii ( 19591 )) of the diffusion equation. This 
solution has been obtained under rather restrictive assumptions that diffusion coefficient 
D{E) and energy loss es b(E) = —dE/dt of the prop agating particles do not depend on time 
t. In our recent work (IBerezinsky fc Gazizovl (120061 ) . to be cited below as paper I) we found 
the analytic solution to the diffusion equation in the expanding universe, valid for the time- 
dependent diffusion coefficient D{E,t) and energy losses b{E,t). The aim of this work is to 
calculate spectra of UHE protons using the BG solution of paper I and to compare them 
with the spectra obtained with the help of the Syrovatsky solution in the above-cited papers. 

The diffusion equation for ultra-relativistic particles propagating in the expanding uni- 
verse from a single source, as obtained in paper I, reads 



n- 



db{E,t) D{E,t) 
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where the coordinate x corresponds to the comoving distance and a{t) is the scaling factor of 
the expanding universe, n = n{t, x, E) is the particle number density per unit energy in an 
expanding volume of the universe, dE/dt = —b{E, t) describes the total energy losses, which 
include adiabatic H{t)E and interaction bint{E,t) energy losses. Qs{E,t) is the generation 
function, that gives the number of particles generated by a single source at coordinate Xg 
per unit energy and unit time. 



According to paper I, the spherically-symmetric solution of Eq. ([T]) is 
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The characteristic trajectory, £' = E'{E, z' 

dE 



is a solution of the differential equation 



— = -[H{t)E + k^t{E,t)]. 



(2) 
(3) 

(4) 
(5) 

(6) 



It gives the energy E' of a particle at epoch z', if this energy is E a.t z = 0; we shall use 
also the notation Eg{E, z) for this quantity. The upper limit Zg in the integral of Eq. ([2]) is 
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provided by maximum energy of acceleration as Eg{E,Zg) = -Emax, or by Zmax whichever is 
smaller. 

To calculate the diffuse flux of UHE protons Jp{E), we sum up contributions of sources 
located in the vertices of a 3D cubic lattice with spacing d. Positions of the sources Xg = 
(^,?7,C) are given by the coordinates ^ = d{i + 1/2), r] = d{j + 1/2), ( = d{k + 1/2), where 
k = 0, ±1, ±2, .. and position of the observer is assumed at ^ = 0, ?7 = 0, C = 0. Thus we 
obtain 

Jp{E) = ^J2''{x^,k,E), (7) 

where n{xijk, E) is given by Eq. ([2]) and 

x.jk = d^/{i + 1/2)2 + + 1/2)2 + (fc + 1/2)2. (8) 



For the propaga ti on of UHE protons in magnetic fields we follow the picture used by 
Aloisio fc BerezinskyI (12004 . 120051 ). namely we assume the turbulent magnetized plasma. 
Magnetic field produced by turbulence is characterized by the value of the coherent mag- 
netic field Be on the basic scale of turbulence Ic, which we shall keep in our estimates as 
1 Mpc. On the smaller scales / < Ic the magnetic field is determined by the turbulence 
spectrum. The critical energy of propagation is determined by the relation r^i^Ec) = Ic, 
where tl is the Larmor radius of a proton. Numerically, Ec = 0.93 x 10^^(-Bc/1 nG) eV. The 
characteristic propagation length in magnetic field is the diffusion length ld{E). It is defined 
as the distance on which a particle is deflected by 1 rad. The diffusion coefficient is defined 
as D{E) = dd{E)/3. For the case ri{E) ^ l^., i.e. when E ^ E^, the diffusion length can 
be straightforwardly found from multiple scattering as 



h{E) 



1.2 



B 



nG 



Mpc, 



(9) 



where Eig = E/{1Q^^ eV) and 5„g = B/{1 nG). M E = E^, h 



At E Ec the diffusion length depends on the spectrum of turbulence. For the Kol- 
mogorov spectrum ld{E) = lc{E / Ec)^^^; for the Bohm regime ld{E) = Ic {E/ Ec). 

The strongest observationa l upper limit on the magnetic fields in our picture is given 
by lBlasi. Buries &: Olintd (119991 ) as -Be ^ 6 nG on the scale Ic = 50 Mpc. In the calculations 
below we assume the representative values of Be in the range (0.1 — 1) nG for Ic = 1 Mpc. 

We do not put as the aim of this paper the d etailed study of diffusion in the time- 
dependent regime. Such a work (lAloisio et al.ll2007bl ) is at present in progress. Here we want 
mainly to demonstrate the validity of the solution for expanding universe obtained in paper I 
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and to perform the numerical comparison of the UHECR diffuse spectra predicted by BG and 
Syrovatsky solutions. The difference is expected to be substantial at energies E < 3x10^^ eV, 
where effects of the universe expansion (in particular, of the CMB temperature growth with 
red-shift) are not negligible. We want also to test the new solution obtained in paper I, 
namely to see the compatibility of the BG and Syrovatsky spectra at high energies, where 
the universe can be considered as the static one, as well as the convergence of BG spectra 
to the universal spectrum, when the source separation d ^ 0. 

The paper is organized as follows: Section [2] addresses the question how reliable is an 
assumption of the diffusion for the low-energy part of UHECR. In Section [3] we calculate the 
diffuse spectra using the time-dependent BG solution. In Section H] we consider the static 
universe, which is the necessary assumption for the Syrovatsky solution, and in Section [S] we 
compare the spectra calculated for the expanding and static universes. The short conclusions 
are presented in Section [61 



2. Why diffusion? 

We argue here that at least at the low-energy end of extragalactic UHECR the diffusion 
propagation is unavoidable for any reasonable magnetic field. We estimate also the magnetic 
field Be for which protons with energies E > 10^'' eV propagate quasi-rectilinearly, hence 
producing the same energy spectrum as in the case of rectilinear propagation. 

To facilitate the calcu lations, let us consider the case of the static universe as in 
Aloisio fc BerezinskyI (120051 ) . namely the universe with 'age' to = Hq^, as follows from 



WMAP observations, and with fictitious 'adiabatic' energy loss of particles b{E) = E Hq, 
where Hq = 72 km/sMpc is the observed Hubble parameter. 

In this picture there is a maximum diffusive distance - magnetic horizon - which is 
determined by the distance traversed by a particle during the age of the universe to: 

rLiE)= r dt D[E,{E,t)] , (10) 
Jo 

where Eg{E,t) is the energy that a particle has at time t, if it is at t = t^. Putting 
dt = -dEg/b{Eg) in Eq. ([TO]), we obtain 



where E^ax = mm[Eg{E, to), E, 



acc 1 
maxj ■ 
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In r^^j.{E) one can recognize (see lAloisio &: BerezinskyI (120051 )) the Syr ovatsky variab le 



X{E,Eg) at Eg = -Emax (for the physical discussion of magnetic horizon see iParizotI (120041 )) 



Let us consider a transition from the diffusive to the rectilinear propagation, allowing 
the considerable deflection angle > 1, when the spectrum is the same as in the rectilinear 
propagation. Two conditions must be fulfilled: 



ld{E) > rhor(^). 



(12) 
(13) 



Eq. ([T^ gives a necessary (but not a sufficient!) condition: the scattering angle on a 
basic scale must be small enough, 6 = lc./ri{E) < 1. It implies the propagation in the regime 
with E > Ec, where ld{E) = 1^ {E/EcY. Considering the low-energy case E <1 x 10^^ eV, 
when the adiabatic energy loss dominates, so that h{E) = EHq and 



-Emax = Eg{E, to) 

one readily obtains from Eq (fTTj) 
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Using Eq. ( 1131) we get 
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From vlIEc) = Ic we obtain Ec = qBclc, where q is an electric charge, which equals 300 for 
a proton, if B is measured in Gauss and E in eV. 



Finally we have 



or, numerically, 



Be < 



E 

qlc 



Br < 1.6 X 10 
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E 



ie' - 1] 
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"1/2 
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-1/2 



nG. 



1017 eV VI Mpc^ 

Therefore, Be < 1.6 x 10^'^ iiG provides the quasi-rectilinear propagation for all protons with 
energies E > 1 x 10^^ eV, while at lower energies the diffusion description is applicable. For 
the reasonably low field Be ~ 0.01 nG the diffusion becomes valid at < 1 x 10^^ eV. 
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3. Diffusive energy spectra of UHECR in the expanding universe 

In the following calculations we shall use the simplified illustrative description of mag- 
netic field evolution with redshift, namely we parametrize the evolution of magnetic config- 
uration (/c, Be) as 

Uz) = IJ{1 + z), B,{z) = 5, (1 + zf-"^, 

where factor (1+2;)^ describes the diminishing of the magnetic field with time due to magnetic 
flux conservation and {l + z)~"^ - due to MHD amplification of the field. The critical energy 
Ec{z) found from riiE) = lc{z) is given by 

EJz)=0.93xlO^^(l + zY-"' 



1 nG 

for /c = 1 Mpc. The maximum redshift used in the calculations is ^max = 4. 

The diffuse flux is calculated for the lattice distribution of the sources (in the coordinate 
space x) with lattice parameter (the source separation) d and a power-law generation function 
for a single source 

where Eq is the normalizing energy, for which we will use 1 x 10^^ eV and go has a physical 
meaning of a source luminosity in protons with energies E > Eq, Lp{> Eq). The correspond- 
ing emissivity Cq = qo/d^, i.e. the energy production rate in particles with E > Eq per unit 
comoving volume, will be used to fit the observed spectrum by the calculated one. 



Using to Eqs. (|2HH]), one obtains the diffuse spectrum as 
c go(7g - 2) 



^tiHq E^ 

^ J [E,{E,z)/Eor^^ exp[-4xfc] dE,{E,z) 

il + z)^nUl + zr + A [4vrA(E,z)]3/2 dE ' ^ > 



where summation goes over the sources like in Eqs. ([7]) and ([8]), the upper limit Zg is provided 
by maximum energy of acceleration as Eg{E, Zg) = E^^^ or by Zm^:^ , whichever is smaller, 7n 



i s the gener ation index and a formul a for dEg/dE can be found in iBerezinsky &: Grigorieva 



( 119881 ) and IBerezinsky et al.l (j2002al ). The analogue of the Syrovatsky variable, \{E,z), is 



given by 

A(E, z) = ^f dz^=^^=== D[Eg{E, z),z]. (21) 
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Fig. 1. — Convergence of the diffusive solution ([7]) to the universal spectrum when the 
distance between sources diminishes from 50 to 10 Mpc shown by numbers on the curves. 
The magnetic field configuration [Be, Ic) = (100 nG, 1 Mpc) with the Kolmogorov spectrum 
of turbulence. The emissivity Cq = q^/d? = 2.4 x 10^^ erg/Mpc^yr is the same for all curves 
The generation spectrum is oc The spectrum of Akeno-AGASA flHonda et a/.lll993 

Shinozakil |2006| ) is shown for comparison. 



First of all we test the correctness of the obtained solution with the help of propa- 
gation theorem (lAloisio &: BerezinskyI |2004| ) . which states that diffusive solution ([7]) con- 
verg es to the univers al spect rum, i.e. one calculated for homogeneous source distribution 
(see iBerezinsky et al.l (l2002al )) when the distances between sources d ^ 0. Fig. [T] demon- 
strates this convergence for the case of strong magnetic field, Be = 100 nG, and the 
Kolmogorov diffusion at E < E^. The diffuse fluxes are calculated using Eq. ([7]) for 
d diminishing from 50 Mpc to 10 Mpc, keeping the same comoving volume emissivity 
-^0 = = 2.4 X 10"^^ erg/Mpc^yr. The generation spectrum is given by 7^ = 2.7. 

One can observe the convergence of the calculated spectra to the universal one when d — > 0, 
in accordance with propagation theorem. 

For small distances between a source and observer, the diffusion approximation for a 
particle propagation is not valid. One can see it from a simple argument that diffusive prop- 
agation time r'^/D must be larger than time of rectilinear propagation, r/c. This condition, 
using D eld, results in r < l^. At distances r < 1^ the rectilinear and diffusive trajectories 
in magnetic fields differ but little and rectilinear propagation is a good approximation as 
far as spectra are concerned. The number densities of particles Q/ATTcr"^ and Q/AttDt, cal- 
culated in rectilinear and diffusive approximations, respectively, are equal at r ~ l^,, where 
Q is the rate of particle production by a source. At smaller distances the rectilinear flux 
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dominates, at larger - the diffusive one. We calculated the number densities of protons 
n{E, r) numerically for both modes of propagations with energy losses of protons taken into 
account, and the distance of transition is taken from equality of the calculated densities. 
We know that this recipe is som ewhat rough and the interp olation between two regimes 
is required, as in calculations by lAloisio fc BerezinskyI (120051 ). However any interpolation 
meets the difficulties by the following reason. 

In fact, the diffusive regime sets up at distances not less than six diffusion lengths 
Id- At distances id ^ ^ < 6ld some intermediate regim e of propag a tion i n magnetic fields 
is valid. When studied in numerical simulations (e.g. Yoshiguchi ( 20031 )). the calculated 
number density n{E, r) satisfies the particle number conservation Airr'^nu = Q, where u is the 
streaming velocity, while in case of interpolation there is only one (unknown) interpolation 
n{ E,r) which conserve s the number of particles. This problem will be studied in detail 
in ( lAloisio et al.l l2007bl ). while for purposes of this paper we can accept the rough recipe 
of transition from diffusive to rectilinear regime as described above. The appearance of 
the artificial pecuharity connected with the accepted propagation transition will be useful 
serving as a mark for the position of transition in the spectrum. 

For rectiline ar propagation for the la ttice distribution of the sources the diffuse spectrum 
is calculated as (IBerezinsky et al.ll2002al ) 
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(22) 



where Cq = qo/d^ = 2.4 x 10"^^ erg/ (Mpc^yr) is the emissivity, Zij^ is the red-shift for a source 
with coordinates i,j, k, and factor (1 + Zijk) takes into account the time dilation. 

The calculated spectra in the expanding universe for Be = I nG and Be = 0.1 nG, 
both for m = 1 and -Emax = 1 x 10^^ eV are show n in Fig. [2] in comparison with all data 
(left panel) and Yakutsk data (lEgorova et al\\2004 ) (right panel). The all -data spectrum is 



obtained using the of energy calibration of all detectors with help of the dip (IBerezinsky et al. 



2002al ). One can observe the peculiarity in predicted spectrum in the right panel [Be = 1 nG) 
at energy E ^2 x 10^^ eV. This is the energy of transition to rectilinear propagation. This 
peculiarity is unphysical and is connected with the simplified description of the transition 
as described above. When magnetic field diminishes a peculiarity shifts toward lower energy 
(left panel) as it should do. At small d the calculated spectra converge to the universal 
spectrum, as it must be. The spectra for very low magnetic field with Be = 0.01 nG is shown 
in Fig. [3l One can notice that the peculiarity is shifted here to the energy E ^ 2x 10^^ eV and 
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Fig. 2. — Diffusive spectra in the expan ding; universe for Be = 0.1 nG (left panel) compared 
witfi Yakutsk data (lEgorova et al\\2004 ) a nd Br. = 1 nG (riglit panel) conapared witfi all-data 



Akenq-AGASA-HiRes- Yak utsk spectrum (IHonda et all Il993l : IShinozakil l2006l : lAbbasi et al. 



2004J : lEgorova et all 12004 ) . The spectra are calculated as a combination of diffusive and 
rectilinear spectra with the Kolmogorov (solid lines) and Bohm (dash lines) regimes of dif- 
fusion at low energies E < Ec. The numbers at the curves indicate the separations of the 
sources in Mpc. The spectral generation index 7^ = 2.7 and maximum acceleration energy 
E^.. = 1 X 10^2 eV. 



above this energy the spectrum is rectilinear. This behavior corresponds to the discussion 
in Section [21 



4. Spectra in the static universe 

In this section we study the Syrovatsky solution of UHECR diffusion equation with aim 
of comparing it with the BG solution for the expanding universe. The Syrovatsky solution 
is valid in case of infinite space with time-independent diffusion coefficient D{E) and energy 
losses b{E). This implies the static universe. 

We define the static universe in which the stationary diffusion equation with the Sy- 
rovatsky solution is embedded in the following way. There is no expansion. The Hubble 
constant is assumed as a formal parameter Hq = 72 km/s Mpc, which defines the "age" 
to of the universe according to WMAP relation Hoto = 0.993, so that cto is the size of the 
universe: the space density of the UHECR sources outside the sphere of radius do is Ug = 0. 
The temperature of CMB photons is constant. To = 2.728 K, and thus the interaction energy 
losses dE/dt = —h{E) are time- independent. 
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Apart from interaction energy losses, we assume the fictitious a diabatic energy losses 
descri bed by dE/dt = —HqE. In this approach we follow the work by lAloisio fc Berezinsky 
fl2005h . 



The universal spectrum in the static universe is different fro m that in the expandin g 



universe. It is defined in the same way as in the expanding universe (IBerezinsky et al. 
namely from the number of particles conservation 



2002aD 



dtQgen[Eg{E,t)] 



(23) 



where Qgcn{Eg) is the generation rate per unit volume and Eg{E,t) is determined by the 
evolution equation dE/dt = —b{E). Note that in the expanding universe Qgcn{E,t) and 
n{E,t) are related to a comoving volume. In contrast to the expanding universe, the ra - 
tio dEg/dE in Eq. fl23p is given by a simple expression (IBerezinsky fc Grigorieval 1 19881 ) 
dEg/dE = b{Eg)/b{E). Using 
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where Cq is the emissivity in particles with energies E > Eq for spectral index 7^ > 2, one 
easily obtains the universal spectrum in analytic form: 
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(25) 
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where -Emax = ^^^[Eg{E ,tQ) , E^^^]. For the diffuse spectra calculations we use the lattice 
distribution of the sources like in the previous section, with the same procedure of transition 
from diffusive to rectilinear propagation, but using for diffusive propagation the Syrovatsky 
solution (see also lAloisio fc BerezinskyI (120051 )). instead of the BG solution. 



In Fig. H] we present the diffuse UHECR spectra calculated, like in Section [3l in the 
grid model with spacing d. UHECR sources are located in vertices of the grid. The spectra 
are calculated as combination of the diffusive flux, described by the Syrovatsky solution, 
combined with the rectilinear flux. They are presented for two magnetic field configurations 
{Bc,lc) equal to (1 nG, 1 Mpc) and (0.1 nG, 1 Mpc) and for different spacings d indicated 
by the numbers in Mpc at the corresponding curves. Like in Section [3], the spectra have 
irregularities at energy of transition from diffusive to rectilinear propagation, caused by the 
rough method of sewing together of the two propagation regimes. 

The calculated diffuse fluxes converge to the universal spectrum (!25l) when d ^ 0, as it 



should be according to propagation theorem. In contrast to calculations by lAloisio fc Berezinsky 
(120051 ). we have obtained the best fit of the data usi ng ■jg = 2.65, which is insignificantly 
different from 7^ = 2.7 of (lAloisio fc BerezinskyI |2005| ). One can see the basic agreement of 
these spectra with those in the expanding universe, including spectrum peculiarities caused 
by transition from diffusive to rectilinear propagation. 
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Fig. 4. — Diffuse spectra in the static universe for Be = 0.1 nG (left panel) and Be = 
1 nG (right panel). The spectra are calculated as combination of diffusive and rectilinear 
spectra with the Kolmogorov (solid lines) and Bohm (dash lines) regimes of diffusion at low 
energies E < E^. The numbers at the curves indicate the separation of the sources in Mpc. 
Observational data are the same as in the right panel of Fig. O 
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5. Comparison of spectra in expanding and static universes 

The direct comparison of the BG and Syrovatsky solutions of diffusion equations is not 
possible because they are embedded in the different cosmological environments. While the 
BG solution is valid for the expanding universe, the Syrovatsky solution is valid only for the 
static universe. Using two different cosmological models for these solutions, there are two 
ways of comparison. 

The first one is given by equal values of parameters in both solutions. In this method 
for BG solution we use the standard cosmological parameters for expanding universe Hq, 
flm, A and maximum red-shift Zmax up to which UHECR sources are still active, magnetic 
field configuration {Bc,lc), separation d and UHECR parameters 7^ and Cq, determined by 
the best fit of the observed spectrum. For the static universe with the Syrovatsky solution 
we use the same parameters Hq, d, {Bc,lc), Ig and Cq. The maximum red-shift in the BG 
solution we chose as that providing the age of universe which equals to to = Hq^ in the 
static universe {zmax = 1.465). This formal method of comparison will be referred to as 
"equal-parameter method". 

Physically better justified comparison is given by the best fit method, in which 7^ and 
£0 are chosen as the best fit parameters for the both solutions, respectively. As a matter of 
fact, we would use the best-fit parameters 7^ and Cq for each solution, if we considered them 
independently. The weakness of this method is a "forced" agreement at the energy range 
of measured UHECR spectrum, provided by the best-fit parameters to the same spectrum. 
Then the difference of the two solutions becomes most appreciable at < 1 x 10^^ eV. 

Equal-parameter comparison of the expanding and static universe solutions are shown 
in Fig. [5] for 7^ = 2.7, £0 = 2.4 x 10^^ erg/Mpc^r and d = 30 Mpc. In the left panel 
i?c = 0.1 nG and in the right panel Be = 1 nG. The universal spectra for the expanding and 
static universes are shown by solid lines. The Kolmogorov diffusion aX E < is presented 
by dash lines and the Bohm diffusion - by dot lines. Note that universal spectra are different 
at < 1 X 10^^ eV as they must be due to excessive energy losses in the expanding universe 
caused by the increasing of CMB temperature T{z) with red-shift. At the highest energy 
end E >Q X 10^^ eV the two solutions coincide exactly because at these energies the energy- 
attenuation time is short, the CMB temperature does not change during time-of-fiight and 
the expanding universe case becomes static. As one can see in Fig. 0, at these energies both 
universal spectra are the same and all spectra for d = 30 Mpc merge into one curve. The 
difference of fluxes at lower energies is naturally explained by the increase of energy losses in 
the expanding universe because of Tcmb{z) dependence, and due to cosmological evolution 
of the magnetic field. 
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Fig. 5. — Equal-parameter comparison of the expanding (exp.) and static (stat.) universe 
solutions for 7^ = 2.7, Co = 2.4 x 10^^ erg/Mpc^yr and d = 30 Mpc. The left panel shows 
the case Be = O.l nG, the right panel is for Be = 1 nG. The universal spectra are shown by 
solid lines. The cases of the Kolmogorov and Bohm diffusion at E < are shown by dash 
and dot lines, respectively. Data in both panels are the same as in the right panel of Fig. [2l 

The best-fit comparison of the expanding and static universe solutions is shown in Fig. El 
The best-fit parameters are (7^ = 2.7, Cq = 2.4 x 10"^^ erg/Mpc^yr) and (7c, = 2.65, £0 = 
5.7 X 10^'* erg/Mpc^yr) for the expanding and static universes, respectively. The left panel 
shows the case Be = 0.1 nG and the right panel Be = 1 nG, both for d = 30 Mpc. One can 
see a good (though "forced") agreement between the expanding universe and static universe 
solutions in the energy range of UHECR observations, which is improved in comparison with 
equal-parameter method because of the choice of the best-fit parameters (7^, Cq) different 
for each case. 

However, we emphasize again that this is a natural way of selection of the solution in 
independent analysis. At energies below 1 x 10^^ eV the larger discrepancies are seen, being 
induced by larger energy losses and evolution of magnetic field in case of expanding universe. 

In conclusion, one can see a reasonably good agreement between the Syrovatsky solution, 
embedded in static universe model, with the BG solution for expanding universe at energies 
E > 1 X 10^^ eV with noticeable discrepancies at smaller energies, which are natural and 
understandable. 
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E, eV E, eV 



Fig. 6. — Best fit comparison of tlie expanding and static universe solutions for d = 30 
Mpc. Best-fit parameters are 7^ = 2.7, Cq = 2.4 x 10^^ erg/Mpc^yr) for expanding and 
7c, = 2.65, Co = 5.7 x 10^^ erg/Mpc^yr) for static solutions, respectively. The universal 
spectrum is shown for the expanding universe. The left panel exposes the case of Be = 0.1 
nG and the right one the case of Be = I nG. The combined data of all arrays are shown for 
comparison. Data in both panels are the same as in the right panel of Fig. [21 



6. Conclusions 



In this paper we study the application of solution of diffusion equation in expanding 
universe obtained in the paper I to the propagation of UHE protons. However, we do 
not consider here the detailed picture with realistic evolution of magnetic field in expanding 
universe, and with a realistic transition from the diffus i ve to t he quasi-rectilinear propagation. 



It is to be considered in our next work lAloisio et al.l (l2007bl ). 



In this paper we demonstrate that solution of diffusion equation found in paper I looks 
quite reasonable when applied to realistic models. Numerically these solutions are similar to 
the Syrovatsky solutions valid for the static universe with the understandable distinctions 
at low energies. 

The diffusion spectra in expanding universe are presented in Figs. [H [2], [3] for the following 
magnetic configurations (BcJc) = (100 nG, 1 Mpc), (1 nG, 1 Mpc), (0.1 nG, 1 Mpc) and 
(0.01 nG, 1 Mpc), respectively. In the latter case the energy spectra at -E > 1 x 10^'^ eV are 
practically the same as in the case of rectilinear propagation. The evolution of magnetic field 
in the expanding universe is given by one illustrative example. The spectra are shown for 
different separations d of the sources. Transition from diffusive to rectilinear propagation is 
described by the simplified recipe, which results in the artificial spectral features as described 
in Section [3l The best fit of the observed spectrum needs the same index of generation 
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spectrum 7^ = 2.7 as in the case of universal spectrum. The diffusive spectrum converges to 
universal spectrum when d ^ 0, as it should be according to the propagation theorem. 

The Syrovatsky solution of the diffusion equation, i.e. one when the diffusion coefficient 
D{E) and energy losses dE/dt = —b{E) do not depend on time t, is valid only for the static 
universe. In the static universe we make several additional assumptions. We introduce 
the Hubble constant Hq as a formal parameter, which determines the fictitious "adiabatic" 
energy losses dE/dt = —EHq. The "age" of the universe to = determines fact the 
sphere of radius tq = cto occupied by the sources. For this universe we calculate the universal 
spectrum given by Eq. fl25l) and the diffusive spectra for given magnetic configuration [Be, Ic) 
and different separations of the sources d. The convergence to the universal spectrum occurs 
when (i — > as it should be. The best fit of the observational da ta is obtained when 7^ = 2.65 , 



which can be considered as a good agreement with the work by lAloisio fc BerezinskyI (j2004j ). 
where 7^ = 2.7 was used. 

For the comparison of the BG and Syrovatsky solutions we use two schemes. In the 
formal one we compare the BG and Syrovatsky solutions for the same magnetic configurations 
(BcJc) and the same emissivities Cq. One can see from Fig. [S] that both solutions coincide 
exactly at E > 6 x 10^^ eV, when effect of universe expansion (most notably variation 
of CMB temperature) can be neglected. At lower energies the difference in the spectra 
naturally emerges due to TcMsit) dependence and thus due to different energy losses in the 
two solutions. Since the energy losses in the expanding universe are larger, the BG spectra 
occur below the Syrovatsky spectra. 

For practical applications the discrepancy in the spectra at energy range 1 x 10^^ —5 x 10^^ 
eV is not essential, because it is eliminated by renormalization of the calculated flux, i.e. 
by changing the emissivity Cq for the static universe solution. In fact, this procedure is 
necessary for fitting of the observed spectra. 

The comparison of two solutions as given above is formal. As was emphasized above, the 
Syrovatsky solution, which is valid for infinite space and time, with time- independent D{E) 
and b{E), needs the specific definition of the static universe. Only in this case one obtains 
the physically viable solution. This solution needs the best fit parameters 7^ and Cq, which 
are different from those in expanding universe. It is physically more meaningful to compare 
the spectra using for static universe its own best fit parameters 7^ and Cq different from that 
in the expanding universe. The comparison shown in Figs. [6] reveals less discrepancies than 
in the formal scheme of comparison. 
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